The Cytokine, Chemokine, and Growth Factor Network of Prenatal Depression

Background: Neuro-immune pathways are engaged in antenatal and postpartum depression. Aims: To determine if immune profiles influence the severity of prenatal depression above and beyond the effects of adverse childhood experiences (ACE), premenstrual syndrome (PMS), and current psychological stressors. Methods: Using the Bio-Plex Pro human cytokine 27-plex test kit, we assayed M1 macrophage, T helper (Th)-1, Th-2, Th-17, growth factor, chemokine, and T cell growth immune profiles as well as indicators of the immune inflammatory response system (IRS) and compensatory immunoregulatory system (CIRS) in 120 pregnant females in the early (<16 weeks) and late (>24 weeks) pregnancy. The Edinburgh Postnatal Depression Scale (EPDS) was used to assess severity of antenatal depression. Results: Cluster analyses showed that the combined effects of ACE, relationship dissatisfaction, unwanted pregnancy, PMS, and upregulated M1, Th-1, Th-2, and IRS immune profiles and the ensuing early depressive symptoms shape a stress-immune-depression phenotypic class. Elevated IL-4, IL-6, IL-8, IL-12p70, IL-15, IL-17, and GM-CSF are the cytokines associated with this phenotypic class. All immune profiles (except CIRS) were significantly associated with the early EPDS score, independent of the effects of psychological variables and PMS. There was a shift in immune profiles from early to late pregnancy, with an increase in the IRS/CIRS ratio. The late EPDS score was predicted by the early EPDS score, adverse experiences, and immune profiles, mainly the Th-2 and Th-17 phenotypes. Conclusions: Activated immune phenotypes contribute to early and late perinatal depressive symptoms above and beyond the effects of psychological stressors and PMS.


Introduction
Antenatal depression (AD) is a worldwide public health problem with negative effects on mother and child [1]. Antenatal depression increases the risk of postnatal depression, low birth weight, operative delivery, preterm birth, preeclampsia, mother and familial

Participants
This research included 120 pregnant women at the Antenatal Care (ANC) Center at Khon Kaen Hospital, a provincial hospital located in Khon Kaen, Thailand. Khon Kaen is a province in Thailand's northeast. Using predetermined inclusion and exclusion criteria, a professional researcher selected the pregnant participants. Pregnant women between the ages of 18 and 49 who resided in the province of Khon Kaen and were able to read and write Thai were eligible. In addition, their pregnancy had not yet reached 16 weeks gestation, and they wanted to get ANC treatment in the hospital until giving birth. Our study is a prospective cohort study whereby we examined two repeated measurements of depression scores and immune data from the early (gestational week < 16) to late (gestational week > 24) prenatal phase. After the first early phase assessment, twentyfour women withdrew, choosing to follow up at alternate facilities. As such, 96 of the 120 women who participated in the early prenatal phase had assessments of depression and immune profiles during the late pregnancy phase. Pregnant women with axis-I mental disorders other than depression-such as bipolar disorder, schizophrenia, autism, substance abuse disorders, psycho-organic disorders, and generalized anxiety disorderwere excluded from the study. We also excluded women with immune, autoimmune, and neuroinflammatory disorders, including inflammatory bowel disease, multiple sclerosis, type 1 diabetes mellitus, systemic lupus erythematosus, rheumatoid arthritis, and psoriasis, as well as women who were taking antidepressants, mood stabilizers, antipsychotic drugs, or immunosuppressive drugs. During the course of the investigation, only supportive psychotherapy, psychoeducation, sleep hygiene education, and family counseling were permitted. Antidepressants, antioxidant supplements, or other medical treatments for depression have not been administered to any of the pregnant women.
Institutional Review Board (IRB) permission was received from both Khon Kaen Hospital and Chulalongkorn University (KEF62036 and COA No. 280/2019). All women provided their informed consent in writing.

Measurements
The sociodemographic data included age, marital status, family income, and the age of the spouse. The women's medical records were examined for their gravidity, parity, pregnancies, body mass index (BMI), and perinatal depression history. The EPDS was used to evaluate depressed symptoms and the degree of depression [28]. The scale consists of 10 items scored from 0 to 3 on a four-point Likert scale, with a total score ranging from 0 to 30. This scale measures the emotional state of pregnant women throughout the previous week. In the Thai version of the EPDS, the validated screening cutoff value for pregnant women was ≥10 (AUC 0.84, sensitivity 60%, and specificity 90%) [29]. This scale has strong reliability and validity, as indicated by research conducted in the upper Northeast area of Thailand among pregnant women in their third trimester [30]. Relationship satisfaction was evaluated using a four-point Likert scale ranging from "satisfied" to "dissatisfied". In addition, we evaluated if the pregnancy was unintended (dummy variable: yes/no). The Abuse Assessment Screen for intimate partner violence (AAS) [31,32] was used to assess partner abuse during pregnancy among our pregnant women. We used item one (AAS1) in the current study since it was the most important in a previous study [7]. The Multidimensional Scale of Perceived Social Support (MSPSS) [33], a 12-item rating measure scored on a seven-point Likert scale (1 to 7 for each item), was used to assess perceived social support [33]. The MSPSS friends subdomain was chosen as the most important domain for this study [7]. The Adverse Childhood Experiences (ACE) Questionnaire was used to evaluate adverse childhood experiences [34]. In the current study, we used [7] psychological abuse, physical abuse, sexual abuse, a formative construct made using the three abuse scores (labeled ACE_abuse), any neglect (mental or physical neglect), a PC extracted from psychological and physical abuse, domestic violence, a family history of substance abuse, and a family history of psychiatric illness (labeled ACE_FDA, ACE-Family-Dysfunction and Abuse [7]. To measure PMS symptoms, the Premenstrual Symptoms Screening Test (PSST) was used. The premenstrual symptoms (items [1][2][3][4][5][6][7][8][9][10][11][12][13][14] and their repercussions in different contexts (A-E: five items) were assessed on a four-point Likert scale ranging from "not at all" to "mild" to "moderate" to "severe" and scored on a scale from 1 to 4 [35]. In the current study, we used the first PC extracted from ten PSST reflecting the affective component of PMS (labelled PSST_depression), the first PC extracted from three physiosomatic PSST symptoms (labelled PSST_physiosom), the first PC extracted from all five PSS interference data (A, B, C, D, and E) (labelled PSST_inference), and the total PSST score [12].

Assays
BD Vacutainer ® EDTA (10 mL) tubes were used to collect fasting blood in the early morning hours. Blood was sampled from 120 women in the early pregnancy phase (<16 gestational weeks) and in 96 of those women in the late pregnancy phase (>24 gestational weeks). We determined the concentrations of cytokines, chemokines, and growth factors in plasma samples using the Bio-Plex Pro human cytokine 27-plex test kit (BioRad, Carlsbad, California, United States of America). The LUMINEX 200 apparatus was used to measure the fluorescence intensities (FI) of the detecting antibodies and streptavidin-PE (BioRad, Carlsbad, California, United States of America). The intra-assay coefficient of variance of all analytes were less than 11 percent. Electronic Supplementary File (ESF), Table S1 lists the names, acronyms, and official gene symbols for all the cytokines/chemokines/growth factors measured in this study. In the current investigation, we chose to perform statistical analyses on the FI data (with the blank analyte removed) since FI values are often a superior alternative to absolute concentrations, particularly when more than one plate is used [21,22]. IL-1β, IL-2, IL-5, IL-7, IL-10, IFN-γ, CCL5, and VEGF exhibited a high number of results below the assay's sensitivity (>60%) and, therefore, were excluded from the analysis performed on the single indicators. Nevertheless, these analytes were used to compute the composite scores when at least 7% (n = 15) of the assays were measurable (thus all, except IL-5). As such, eight analytes were excluded from the single indicator analyses, but IL-1β, IL-2, IL-7, IL-10, IFN-γ, CCL5, and VEGF could be used as prevalences in the computation of composite scores. ESF, Table S2 provides a summary of the immunological profiles studied in this investigation [20][21][22]. Where needed, cytokines were used after logarithmic (log10), square root, fractional rank-based normal transformations, or Winsorization, and all data were processed as z transformations.

Data Analysis
We used analysis of variance (ANOVA) to compare scale variables across research groups and Chi-square (χ 2 ) or Fisher's exact probability test to compare categories. Using Pearson's product-moment correlation coefficients and the point-biserial correlation coefficient, correlations between scale variables or scale and binary variables were analyzed. Using Generalized Linear Model (GLM) analysis, the connection between the EPDS score and immune profiles was investigated. An automatic stepwise binary logistic regression analysis was performed to determine the most accurate predictors of prenatal depression versus the absence of antenatal depression using the immune profiles, PSST and ACE scores, relationship satisfaction, planned versus unplanned pregnancy, AAS1, and the MSPSS_friends score. Using multiple regression analysis, the most influential factors (immune and psychological variables and PMS) affecting the EPDS score were determined. In addition to the manual regression technique, we used an automated technique with p-values of 0.05 for model entrance and 0.10 for model exclusion. We determined the model statistics (F, df, and p values) and the total variance explained (R 2 ), as well as the standardized β coefficients with t statistics and p-values for each predictor. Additionally, the variance inflation factor and tolerance were evaluated to identify any collinearity or multicollinearity concerns. Heteroskedasticity was determined with the use of the White and modified Breusch-Pagan homoscedasticity tests. Generalized estimating equations (GEE), repeated measures, were employed to examine the associations between the repeated assessments of the EPDS and biomarker data. For statistical significance, two-tailed tests with a p-value of 0.05 were utilized. The various immune profiles were investigated using principal component analysis (PCA) to check whether one general immune PC could be extracted. The KMO statistic (cut-off value > 0.6), Bartlett's test, and the anti-image correlation matrix were used to evaluate the data's factoriability. To be accepted as a valid PC, the variance explained should be >50.0% and all loadings should be >0.7 on this first PC. We used two-step cluster analysis to delineate a cluster of patients with increased EPDS scores and immune profiles, both entered as continuous variables, and the antenatal depression diagnosis entered as a categorical variable. The number of clusters was determined automatically, and the cluster solution was evaluated using the silhouette measure of cohesion and separation (threshold value > 0.5). All statistical analyses were conducted using version 28 of IBM, SPSS for Windows.

Results
Results of cluster analyses performed on the early phase data. We found that the early phase EPDS score was strongly associated with the total PSST score (r = 0.688, p < 0.001, n = 120), relation dissatisfaction (r = 0.406, p < 0.001), unplanned pregnancy (r = 0.413, p < 0.001), and AAS1 (r = 0.403, p < 0.001). A large part (31.0%; F = 13.34, df = 3/111, p < 0001) of the variance in the total EPDS score was explained by the combined effects of unplanned pregnancy (p < 0.001), relation dissatisfaction (p = 0.017), and AAS1 (p < 0.001). Since these psychological factors explain such a large part of the variance, we also computed residualized EPDS (resEPDS) values after covarying for relation dissatisfaction, unplanned pregnancy, and AAS1. This new EPDS score thus reflects severity of depression independent of the effects of adverse psychological stressors.
A two-step cluster analysis was performed with the diagnosis of antenatal depression (EDPS ≥ 10) as categorical factor, total EPDS and resEPDS scores, and the three major cytokine profiles (M1, Th-1, and Th-17) as continuous variables. A two-cluster solution yielded the best solution with a silhouette measure of cohesion and separation of 0.68, whereby the first cluster comprised 76 people and the second cluster comprised 44 people. Figure 1 shows that both EPDS scores and all immune profiles, except CIRS, were significantly higher in cluster two than in cluster one. ESF, Table S3 shows the results of univariate GLM analyses of the immune profiles with age and BMI as covariates. The top three most important profiles of cluster two were (in descending order): Th-1, IRS, and T cell growth. ESF, Table S4 shows the outcome of univariate GLM analyses performed on the single cytokines-growth factors. IL4, IL-6, IL-8, IL-12 p70 , IL-15, IL-17, and GM-CSF were significantly higher in cluster two than in cluster one. ESF, Figure S1 shows a clustered bar graph with the differences in those cytokines.  Table 1 shows the socio-demographic and clinical data in both clusters. The range of the early EPDS score was 0-20, and the range of the late EPDS score was 0-16. There were no significant differences in age, pre-pregnancy BMI, gravidity, parity, and a history of Caesarian delivery between both clusters. Unplanned pregnancy, relation dissatisfaction, and AAS1 were all significantly higher in cluster two than in cluster one. PSST_depression, PSST_physiosom, PSST_interference, and PSST_total were significantly higher in cluster two than in cluster one, whilst MSPSS_friends was significantly lower in cluster two. Psychological abuse, physical abuse, total abuse, ACE1, domestic violence, and family dysfunction were significantly higher in cluster two than in cluster one. Binary logistic regression analysis showed that PSST_interference (Wald = 6.02, df = 1, p = 0.014), PSST_depression (Wald = 5.30, p = 0.005), and All_abuse (Wald = 9.17, p = 0.002) predict cluster two versus cluster one (χ2 = 43.57, df = 3, p < 0.001) with a Nagelkerke value of 0.416. The total number of correctly classified patients was 84.2% (sensitivity = 77.3% and specificity = 88.2%). As such, the cluster analysis performed on the early pregnancy data has generated a cluster of pregnant women with increased adverse experiences, immune, and depression scores.
ESF, Table S5 shows the results of GEE (repeated measurements) analyses with the early and late EPDS and immune profiles as dependent variables, and time (early-late pregnancy phase), the early pregnancy clusters, and their interaction as explanatory variables. There was a significant effect of the interaction pattern time × diagnosis (and significant time and diagnosis effects) on the repeated EPDS assessments. We found a significant increase in cluster one people (p = 0.010), but in cluster two there was a significant decrease in the EPDS score from the early to the late phase (p < 0.001). The same table shows the increased late EPDS values in cluster one are accompanied by increased Th-1, Th-17, IRS, growth factor, and T cell growth profiles. The lowered late EPDS scores in cluster two are accompanied by lowered Th-1 and Th-17 profiles.

Immune Profiles of Early Prenatal Depression
The first major findings of this study are that (a) all immune profiles (except CIRS) during early pregnancy were significantly linked with the EPDS score, independent of the effects of psychological variables and PMS; and (b) elevated CXCL8 (IL-8), IL-15, IL-17 and GM-CSF are the primary cytokines associated with prenatal depression, while IL-4, IL-6 and IL-12 associate with an increased EPDS score in cluster 2 participants. Previous research has shown that prenatal depression is accompanied by elevated C-reactive protein, a positive acute phase protein, and decreased zinc, a negative acute phase reactant, suggesting a modest inflammatory response [1,18]. Our findings that IL-6 is linked with the severity of prenatal depression are consistent with those of previous studies [19,36,37]. These authors also found that IL-1β [19] and TNF-α [34] may be related with moderate depression symptoms in early pregnancy to midgestation. In contrast to the findings of Sha et al. [19], we discovered that CXCL8 (IL-8) was associated with the severity of the condition. Increased CXCL8 was also detected by Corwin et al. [38] who reported that postpartum depression might be predicted by a high IL-8/IL-10 ratio. Karlsson et al. [39] found that although IL-6 and TNF-α were unrelated to prenatal symptoms of depression, IL-12, IL-9 and IL-13, as well as the IFN-γ/IL-4 ratio and IL-5, were all positively associated with depressive symptoms.
In addition, Maes et al. [40] revealed that pregnant women with a history of major depression had elevated levels of IL-6 and soluble IL-1 receptor antagonist (sIL-1RA), suggesting sensitization of the immune-inflammatory response. In accordance, Christian et al. [41] observed that blood levels of macrophage migration inhibitory factor (MIF) were higher after one week of influenza vaccination in pregnant women with greater depressive ratings, suggesting immunological sensitization. Research has connected prenatal depression to inflammatory pregnancy morbidities such as preeclampsia, preterm birth, and gestational diabetes [42][43][44][45][46].
Current hypotheses and reviews demonstrate that immunological pathways that contribute to non-pregnancy major depression also contribute to perinatal depression. Hence, Leff-Gelman et al. [47] postulated that the same immune-inflammatory pathways linked with major depression also operate in the placenta during prenatal depression. Activation of Toll-Like Receptors (TLR) in placental immune cells, for instance, may promote a Th1/Th17 shift with an increase in pro-inflammatory cytokines that cause depressive symptoms [47]. Our analysis demonstrates, however, that prenatal depression is accompanied by hyperactive Th-1 as well as Th-17 phenotypes. The comprehensive review of Roomruangwong et al. [1] reported that not all neuro-immune and neuro-oxidative pathways that define non-pregnancy-related depression have a role in perinatal depression, including increased bacterial translocation. The pregnancy-specific effects of sex hormones (oestrogen and progesterone) and the significant activation of the hypothalamic-pituitary-adrenal (HPA) axis during pregnancy may explain some of the discrepancies in the literature, as explained below.

IRS and CIRS Activity throughout the First Trimester
The second major finding is that, as with non-pregnancy-related MDD [20][21][22], prenatal depression severity is linked with greater IRS vs. CIRS responses, despite the presence of a robust Th-2 response and elevated IL-4 levels. The latter may impose negative feedback on M1 and Th-1 phenotypes [20]. Moreover, elevated glucocorticoid levels during pregnancy may attenuate the synthesis and release of Th-1 and M1 cytokines, encouraging a Th-2 shift [48]. Thus, the observed levels of IRS (M1, Th-1, and Th-17) and CIRS (Th-2 and Treg) cytokines throughout various phases of pregnancy are dictated by the tight balance between IRS and CIRS activity and elevated glucocorticoid levels. In the past, we have hypothesized that enhanced HPA-axis activity during pregnancy may reduce IRS activation during pregnancy, while boosting Th-2 and Treg activity [1]. Complicating matters, there is a negative correlation between inflammatory markers and cortisone levels in healthy pregnant women, but not in those with prenatal depression [49,50], whereas other authors hypothesized that dysfunctions of the HPA-axis may be mediated by immunological and placental mechanisms [51]. Pregnancy-specific alterations in sex hormones, particularly progesterone, may further modulate the homeostatic setpoint between the IRS and CIRS phenotypes [1].
These hormonal changes may explain why the IRS pathway alterations seen in this study are substantially less profound than those observed in non-pregnant MDD/BD patients [1,20,22]. These hormonal effects may also explain why not all research on proinflammatory characteristics in prenatal or perinatal depression had positive outcomes. For example, a first investigation was unable to find a link between IL-6, soluble IL-6 receptor (sIL-6R), and sgp130 (thus measuring IL-6 signaling, including trans-signaling), sIL-1RA, and leukaemia inhibitory factor and postpartum depression, despite substantial correlations with postpartum blues [52]. A pilot study reported an inverse connection between IL-1β, TNF-α, and IL-17 and perinatal depression [53]. TRAIL (TNF-α-related apoptosis inducing ligand), M-CSF (macrophage colony-stimulating factor), and CX3CL1 (fractalkine) were reduced in prenatal depression, according to Edvinsson et al. [54].

Variations in Immune Profiles between Early and Late Pregnancy
The third major finding of this study is that there was a shift in immune profiles from early to late pregnancy, characterized by a sharp increase in the IRS/CIRS ratio due to a decrease in CIRS and increases in the M1 and growth factor profiles, indicating that the other profiles (Th-1, Th-2, Th-17, and CIRS) remained more stable over time. Earlier, Blackmore et al. [55] demonstrated that cytokines, particularly IL-6 and TNF-α, are remarkably stable in pregnant women. However, the actual changes in EPDS from early to late pregnancy were positively associated with Th-2 and Th-17 phenotypes, indicating that again, the balance between different types of IRS and CIRS (as well as glucocorticoids and sex hormones) impact the severity of late prenatal depression, and consequently postnatal depression, which is strongly predicted by prenatal depression [1].
Pregnancy is known to cause alterations in the immune system's equilibrium. The invasion of the uterine epithelium by the blastocyst during implantation triggers an immuneinflammatory response [56][57][58], in which endometrial cells and immune cells alike produce increased quantities of IL-6, TNF-α, and CXCL8 [59]. In the second phase, the immune system gains a more anti-inflammatory mode [58], as seen by enhanced CIRS responses in our research. Induction of an inflammatory response in the latter days of pregnancy facilitates uterine contractions, birth, placental rejection, and protection against puerperal sepsis [58,60]. Hence, symptoms of enhanced peripheral inflammation are present towards the end of a full-term pregnancy but are kept in check by heightened immunoregulatory feedback systems, which are achieved at least in part by higher HPA-axis activity and progesterone levels, as well as immunological homeostatic mechanisms [1]. As such, the immune pathophysiology of prenatal depression is not comparable with that of nonpregnancy-associated depression.

Results of Cluster Analysis
The fourth major finding of this study is that we were able to identify a cluster of early pregnancy women with elevated EPDS scores, MCAS/PMS, immunological profiles, adverse childhood experiences, and elevated current psychological stressors, such as unexpected pregnancy and intimate partner dissatisfaction. As a result, these early and current adverse life experiences, PMS, depressive symptoms, and IRS activation shape a cluster of pregnant women that is well-defined. Notably, the latter class is statistically more significant than the diagnosis of prenatal depression as determined by an EPDS ≥ 10 cutoff score. Moreover, this new phenotypic class was also characterized by increased levels of pro-inflammatory (IL-6, IL-8, IL-12, IL-15, IL-17, and GM-CSF) and anti-inflammatory (IL-4) cytokines. Our results suggest that the combined effects of these psychological (ACE and acute psychological stressors) and biological (immune profiles and PMS) factors, plus the ensuing early pregnancy-related depressive symptoms, form a new phenotypic class. The latter combines a lifetime trajectory of psychological (ACE) and hormonal stressors (PMS) and neuro-immune pathways (PMS) with the ensuing prenatal depressive symptoms and activated immune networks.
In the current investigation, however, we determined that these factors yielded cumulative effects and that there were no significant connections between ACE or PMS and the immunological profiles of prenatal depression. In contrast to earlier research, which indicated that ACE, PMS or menstrual distress, and current psychological stressors impact indicators of immunological IRS activation, the current study revealed no such association. First, ACE were shown to be highly related with the stimulated production of cytokines, chemokines, and growth factors by stimulated whole blood culture supernatant, but not unstimulated whole blood culture supernatant, suggesting immunological sensitization to immune injury [21,22]. In the current investigation, however, we examined plasma levels of the immunological products, which are more closely related to unstimulated levels in culture supernatant than stimulated levels. Second, PMS and menstrual distress are related to elevated levels of proinflammatory chemokines, acute phase reactants, and tryptophan catabolites, which are triggered by inflammatory processes [13,25,61]. Third, ACEs and psychosocial stresses in adult life are known to stimulate immune-inflammatory cytokines including IL-6 [21,26,27]. Even modest psychological stressors, such as examination stress, may stimulate the production of pro-inflammatory cytokines, including IL-6, TNF-α and IFN-γ [62]. Nevertheless, such changes in the prenatal period are likely obscured by the impact of the greatly elevated HPA-axis activity and progesterone levels, which alter the homeostatic equilibrium between IRS and CIRS phenotypes from early to late pregnancy [1]. Our results that psychological stressors and PMS had less of an influence on the late EPDS score than the early EPDS score are similarly suggestive of such confounding effects.

Limitations
This work would have been more intriguing if we had assessed HPA-axis activity and plasma progesterone levels, as well as cytokines, chemokines, and growth factors, in the supernatant of stimulated whole blood cultures.

Conclusions
Independent of the impact of psychosocial factors and PMS, all immunological profiles, including M1, Th-1, Th-2, and IRS, were significantly correlated with the early EPDS score. Cluster analysis revealed that the combined impacts of ACE, relationship dissatisfaction, undesired pregnancy, PMS, and activated immunological profiles (including elevated levels of IL-4, IL-6, CXCL8, IL-12, IL-15, IL-17, and GM-CSF), as well as the resulting early depressive symptoms, define a new endophenotypic class. Immune profiles changed from early to late pregnancy, with a significant rise in the IRS/CIRS ratio. Immune profiles, specifically the Th-2 and Th-17 phenotypes, predicted the actual changes in EPDS score from early to late pregnancy. Activation of the immune system contributes to depressive symptoms during pregnancy.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/ 10.3390/brainsci13050727/s1, ESF, Table S1. Cytokines, chemokines and growth factors measured in the current study; ESF, Table S2. Computation of the immune profile composite scores assessed in this study; ESF, Table S3. Results of univariate GLM analyses of the immune profiles and the two-step-cluster analysis-generated classes with age and BMI as covariates; ESF, Table S4. Results of univariate GLM analyses of the single cytokines and growth factors and the two-step-cluster (TSC)-analysis-generated classes with age and BMI as covariates; ESF, Figure S1. Results of cluster analysis performed on the Edinburgh Postnatal Depression Scale scores and immune profiles; ESF, Table S5. Results of general estimation equations (repeated measurements) analyses with the early and late Edinburgh Postnatal Depression Scale (EPDS) scores and immune profiles as dependent variables and time (early-late pregnancy phase), the cluster analysis-generated classes, and their interaction as explanatory variables.

Data Availability Statement:
The dataset generated during and/or analyzed during the current study will be available from the corresponding author (MM) upon reasonable request and once the dataset has been fully exploited by the authors.